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This report summarizes our efforts on the funded project “Quality Control and Analysis 
of Microphysical Data Collected in TRMM Airborne Validation Experiments”, NASA 
NAG5-9663, Andrew Heymsfield, P. I. We begin this report by summarizing our 
activities in FY2000-FY2004. We then present some highlights of our work. The last part 
of the report lists the publications that have resulted from our funding through this grant. 

SUMMARY OF ACTIVITIES 

In FY2000-2001, we focused on analyzing data from the TRMM field campaigns 
TEFLUN-B, LBA and Kwajex. Seminars reporting on the results of our analysis were 
presented by Andy Heymsfield at NASA GSFC, Colorado State University, and at the 
National Center for Atmospheric Research. A talk describing the results was presented at 
the European Geophysical Society (EGS) meeting in Nice. 

In FY2001-2002 , progress was made in five areas. We completed a manuscript that 
described the microphysical properties of the stratiform ice cloud regions observed 
during the TRMM field campaigns. (See Reference List, Heymsfield et al., 2002). The 
study also developed parameterizations for these properties for use by the remote sensing 
and modeling communities, as described in the section on “Primary Findings” below. The 
results were presented at the TRMM Workshop in October 2001. We also worked with 
Jeff Stith to incorporate some of our findings in a paper, which he published (Stith et al., 
2002 ). 

The data from the stratiform ice cloud regions were used to study the relationship 
between the radar reflectivity at various radar wavelengths and the precipitation rate, for 
ice clouds. The wavelengths include those of the TRMM radar and those proposed for the 
Global Precipitation Mission (GPM). The results of this work were presented at the 
TRMM Science Meeting during July 2002. 

The data from the stratiform ice cloud regions were used to characterize the fall velocity 
characteristics of tropical ice cloud populations. Two papers describing these activities 
appear in the list of references (Heymsfield, 2003). 


We assisted with the software development effort for the Common Microphysics Product 
(CMP) definition. Specifically, we helped to debug the software that was developed by 
Alexei Korolev, and assisted Jeffrey Stith and Julie Haggerty with the quality control. 
This work led to a publication (Kingsmill et al., 2004, see below). 

In FY2003-2004, we continued our work on a project that related properties of the 
particle size distributions and particle habits during the periods defined by the CMP effort 
to coincident measured radar reflectivity. We derived gamma fit coefficients that 
described the intercept, slope, and dispersion that characterize the particle size 
distributions. These size distributions were related to those measured in hurricanes during 
CAMEX-4 and also in midlatitude storms. A publication describing this work appears in 
the reference list (Heymsfield et al., 2004). 

We also worked on a characterization of the properties of the particle size distributions in 
convective cells sampled by the Citation aircraft during LB A and KWAJEX. This and 
associated work led to a publication (Stith et al., 2004). Seminars describing our activities 
were presented at NCAR, Colorado State University, and at Cal Tech University. 

PRIMARY FINDINGS 

Some of the main findings from our analysis of the TRMM microphysical data were 
reported in the articles identified in the reference list. The measurements of particle size 
distributions (PSDs) and ice particle habits we helped to acquire during slow, 
Lagrangian-type spiral descents through deep subtropical and tropical stratiform cloud 
layers in Florida, Brazil and Kwajalein, Marshall Islands were examined and synthesized. 
Our purpose with this type of flight pattern was to learn more about how the PSDs 
evolved with height in the cloud layer, and to develop parameterizations for the resulting 
observations. Most of these spirals began at cloud top, with temperatures (T) as low as - 
50C, and they ended at cloud base or below the melting layer (ML). Our analysis has 
shown that the PSDs broadened downwards, with the diameters of the largest particles 
increasing from several millimeters at cloud top to one centimeter or larger at the top of 
the ML. There was some continued particle growth in the ML. The concentrations of the 
sub-1 mm size particles decreased with vertical distance from the top to bottom of the 
cloud layer. The result was a consistent change in the PSDs in the vertical. Aggregation, 
as ascertained from both the changes in the PSDs and evolution of particle habits, was 
responsible for these trends. 

Until recently, we found that the PSDs were generally well represented by exponential 
curves. However, our conclusions are changing, and the resulting analysis is taking quite 
a bit of time. We searched for a better representation of the PSDs over the size range 
from tens of microns to few cm than could be provided by exponentials, as the latter did a 
poor job of representing the PSDs below 500 or so microns. We are progressing very well 
on this issue, and are now using a gamma distribution with a variable order that depends 
primarily on temperature. 
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The use of exponential size distributions did reveal some very interesting results. The 
slope X of the exponentials varied systematically with temperature, and conformed to 
previous studies for mid-latitude frontal and cirrus layers. However, the temperature did 
not appear to be the fundamental factor controlling X. Rather it depended on such factors 
as the collection kernel for aggregation and the time available for this process to occur. A 
direct relationship was found between X and the intercept No .This relationship was not 
found clearly in earlier studies because the PSDs in our observations had better fidelity 
with less scatter. In the ice regions and upper part of the ML, the N 0 and X decreased, 
and in the lower part of the ML they began to increase as the largest particles melted. 
Figure 1 shows an example of the N 0 vs. X relationship assuming an exponential 
distribution; this will change somewhat when our investigation of the gamma functional 
form is completed. 
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Figure 1 N 0 vs. X relationship for spirals on six days during KWAJEX. Different colors 
represent different spirals. Large symbols represent data points at temperatures of 0 and 
3°C 


General expressions were developed relating N 0 to X for various bulk 
microphysical and radar variables, as shown below in Table 1 . These expressions are 
undergoing minor revision to take into account our observations that the PSDs were 
better represented by gamma distributions than by exponentials (gamma of order 0). 
Nevertheless, we feel that it is useful to show the types of relationships that we 
developed. Relationships are included in Table 1 to obtain ice water content (IWC), radar 
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reflectivity factor (Z), mass and reflectivity-weighted fall velocities (Vm and Vz), 
precipitation rate (R), cross sectional area of the particle population (A), and effective 
radius (r e ). These expressions appear to apply to both tropical and mid-latitude clouds 
and facilitate the specification of N 0 and k for cloud and climate models. These 
expressions are also useful for interpreting vertically pointing radar data. Given N 0 and k 
other properties, e.g., mass-weighted fallspeed, can be derived. The approach is well 
suited to calculation of IWC, precipitation rate, and radar-related parameters where 
reflectivities are from 0 to 25 dBZ e . 
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Tkble 1 

Equations to Derive N 0 , A from Microphysical and UmAm* Variables 
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* Incomplete gamma function: 7 ( 0 , 1 ) = *)e l dt. Units milwm otherwise stated: cgs. 

Equations based on the following relationships: 

Mass: m = f p e D s . Effective density: p e = l:(A r ) n I> a > k = 0.4, n = 1.5, a = -0.5 
Area Ratio: Ar = aV*, a = 0.40, 6 = -0.09 
Terminal Velocity: Vt = ~ “ ^Ar" ~ ^ 6/ . 

Drag Coefficient Parameters: aj — 1.653, bf — 0.465. Kinematic viscosity: 1 /. 

Air density: p fl 
K - fc/[6(« - 1) + 3 + <*] 

C =a > /<«- 1 )a / (^)*/ v ( 1 - 2 ^), P = 500 hPa, T = 0*C 


The data used in this analysis was incorporated into calculations reported in the two 
papers on terminal velocity by Heymsfield (2003). 
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